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Q , Abstract 

f \ \ In this paper we present a new approach for studying the dynamics of spatially inhomogeneous 

cosmological models with one spatial degree of freedom. By introducing suitable scale-invariant 
dependent variables we write the evolution equations of the Einstein field equations as a system of 

rs, autonomous partial differential equations in first-order symmetric hyperbolic format, whose explicit 

^ ■ form depends on the choice of gauge. As a first application, we show that the asymptotic behaviour 

near the cosmological initial singularity can be given a simple geometrical description in terms of 
the local past attractor on the boundary of the scale-invariant dynamical state space. The analysis 
suggests the name "asymptotic silence" to describe the evolution of the gravitational field near the 
cosmological initial singularity. 
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1 Introduction 

The simplest cosmological models are the Friedmann-Lemaitre (FL) cosmologies, which describe an 
expanding Universe that is exactly spatially homogeneous and spatially isotropic. It is widely believed 
that on a sufficiently large spatial scale the Universe can be described by such a model, at least since the 



5-j \ time of last scattering of primordial photons with unbound electrons. 

There are, however, compelling reasons for studying cosmological models more general than the FL 
models. Firstly, the observable part of the Universe is not exactly spatially homogeneous and isotropic on 
any spatial scale and so, from a practical point of view, one is interested in models that are "close to FL" 
in some appropriate dynamical sense. The usual way to study deviations from an FL model is to apply 
linear perturbation theory. However, it is not known how reliable the linear theory is and, moreover, 
in using it one is a priori excluding the possibility of finding important non-linear effects. Secondly, it 
is necessary to consider more general models in order to investigate the constraints that observations 
impose on the geometry of spacetime. Thirdly, it is important to classify all possible asymptotic states 
near the cosmological initial singularity (i.e., near the Planck time) that are permitted by the Einstein 
field equations (EFE), with a view to explaining how the real Universe may have evolved. 

For these and other reasons it is of interest to consider a spacetime symmetry-based hierarchy of 
cosmological models that are more general than FL. On the first level above the FL models are the spatially 
homogeneous (SH) models, i.e., models which admit a 3-parameter isometry group acting transitively on 
spacelike 3-surfaces, and expand anisotropically. This class has been studied extensively, and a detailed 
account of the results obtained up to 1997 is contained in the book edited by Wainwright and Ellis (WE) 
| p3[ . On the second level of the hierarchy are cosmological models with two commuting Killing vector 
fields (i.e., models which admit a 2-parameter Abelian isometry group acting transitively on spacelike 2- 
surfaces), which thus admit one degree of freedom as regards spatial inhomogeneity. This class of models, 
which are referred to briefly as Gi cosmologies , are the focus of the present paper. In generalising from 
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1 INTRODUCTION 2 

SH cosmologies to G2 cosmologies one makes the transition from ordinary differential equations (ODE) to 
partial differential equations (PDE) in two independent variables as regards the evolution system of the 
EFE, with the inevitable increase in mathematical difficulty. For both classes of models one has available 
the four standard methods of systematic investigation: 

(i) derivation and analysis of exact solutions, 

(ii) approximation methods of a heuristic nature, 

(iii) numerical simulations and experiments, and 

(iv) rigorous qualitative analysis. 

All four methods have been used to study Gi cosmologies with varying degrees of success, subject to 
significant limitations. 

We now give a brief history of Gi cosmologies. To the best of our knowledge, the first development 
was the study by Gowdy of a class of solutions of the vacuum EFE with compact space sections and 
an Abelian Gi isometry group, now called Gowdy spacetimes p4, E5[. Although vacuum, they can be 
regarded as idealised cosmological models because they have a preferred timelike congruence, start at 
a curvature singularity, and either expand indefinitely or recollapse. These solutions could represent 
the early stages of the Universe during which the energy-momentum-stress content is not dynamically 
significant. As regards Gi cosmologies with a perfect fluid matter source, the earliest paper was by Liang 
p5| , who used approximation methods to study the evolution of matter density fluctuations. A variety of 
exact perfect fluid G2 cosmologies have been discovered, starting with Wainwright and Goode |54), and 
more recently by Ruiz and Senovilla |4J| , Mars and Wolf [J58| and Senovilla and Vera |4g, [48[ . Most exact 
solutions have been derived by imposing a separability assumption on the metric components, so that 
the EFE decouple into two sets of ODE. As regards numerical simulations, work began in the 1970s (see, 
e.g., Centrella and Matzner [Q). Recent work (e.g., Berger and Moncrief [g| and Berger and Garfinklc 
J7J) has focussed on investigating the nature of the cosmological initial singularity in Gowdy vacuum 
spacetimes. Rigorous qualitative analysis has also focused on the past asymptotic behaviour of these 
spacetimes, starting with the paper of Isenberg and Moncrief |3^] on the diagonal subcase. Recent work 
by Kichenassamy and Rendall |J4J and by Anguige g has considered the Gowdy vacuum spacetimes 
with spatial topology T 3 and diagonal perfect fluid Gi cosmologies, respectively. We also refer to Rein 
[ pd[ for related results for a different matter model. 

In summary, almost all of the research using methods (iii) and (iv) above has focused on the Gowdy 
vacuum spacetimes. The mathematical reasons for this choice are twofold: vacuum Gi models are much 
more tractable than non- vacuum ones, and the assumption of compact space sections makes numerical 
simulations easier since it avoids the problem of boundary conditions at spatial infinity. These works are 
nevertheless of considerable physical interest in view of a conjecture by Belinskii, Khalatnikov and Lifshitz 
(BKL) that "matter does not matter" close to the cosmological initial singularity, i.e., that matter is not 
dynamically significant in that epoch (see Lifshitz and Khalatnikov |36 , p200, and Belinskii et al H, 
p532 and p538). We shall refer to this conjecture as BKL I. 

In this paper we focus on G2 cosmologies with perfect fluid matter content, incorporating vacuum 
models as an important special case. The overall goal is to provide a flexible framework for analysing the 
evolution of these models in a dynamical systems context. Our approach, which uses the orthonormal 
frame formalism ,[J has three distinctive features: 

(i) first-order autonomous equation systems, 

(ii) scale-invariant dependent variables, 

(iii) evolution equations that form a system of symmetric hyperbolic PDE. 

Dynamical formulations employing (i) and (ii) have proved effective in studying SH cosmologies (see WE 
|J53| , Ch. 5, for motivation). We expect similar advantages to be gained in the study of G2 cosmologies. 
In the SH case the scale-invariant dependent variables are defined by normalisation with the volume 

x See, e.g., MacCallum |37[, and for an extended set of equations, van Elst and Uggla pa]. 
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expansion rate of the G3-orbits, i.e., the Hubble scalar H . In the present case, however, we define scale- 
invariant dependent variables by normalisation with the area expansion rate of the G2-orbits.1l in order 
to obtain the evolution equations as a system of PDE in first-order symmetric hyperbolic (FOSH) format. 
In this way we make available an additional set of powerful analytical tools, that ensures local existence, 
uniqueness and stability of solutions to the Cauchy initial value problem for Gi cosmologies and provides 
methods for estimating asymptotic decay rates.n The FOSH format also provides a natural framework 
for formulating a concept of geometrical information propagation, by which we mean the propagation at 
finite speeds of jump discontinuities in the initial data set.n 

We now digress briefly to describe some aspects of SH dynamics. The use of scale- invariant dependent 
variables to study SH models led to an important discovery, namely that self-similar solutions of the 
EFE play a key role in describing the dynamics of SH models, in that they can approximate the early, 
intermediate and late time behaviour of more general models. We refer to WE |Q, Ch. 5, for details 
and other references. A self-similar solution admits a homothetic vector field, which, in physical terms, 
means that as the cosmological model expands, its physical state differs only by an overall change in the 
length scale, i.e., the dynamical properties of the model are scale-invariant. The above discovery had 
been anticipated some years earlier by Eardley |13) , who observed that SH models of Bianchi Type—I, 
while not self-similar, are asymptotically self-similar. By this one means that in the asymptotic regimes, 
i.e., near the cosmological initial singularity and at late times, their evolution is approximated by self- 
similar models. In other words, these simple models have well-defined asymptotic regimes that are scale- 
invariant. In general, however, SH cosmologies are not asymptotically self-similar. For example, the 
well-known Mixmaster models (vacuum solutions of Bianchi Type-IX; see Ref. |4(|) oscillate indefinitely 
as the cosmological initial singularity is approached into the past, and thus do not have a well-defined 
asymptotic state (see e.g. Ref. H and WE p3|, Ch. 11). Nevertheless, as follows from the dynamical 
systems analysis, the Mixmaster models are successively approximated by an infinite sequence of self- 
similar models (Kasner vacuum solutions) as they evolve into the past towards the cosmological initial 
singularity. The mathematical reason for the above phenomena is that the self-similar solutions arise as 
equilibrium points (i.e., fixed points) of the evolution equations. These equilibrium points, in conjunction 
with the Bianchi classification of the G3 isometry group, determine various invariant submanifolds of 
increasing generality that provide a hierarchical structure for the SH dynamical state space. In other 
words, the self-similar solutions play a key role as building blocks in determining the structure of the 
SH dynamical state space. We anticipate that self-similar models will play an analogous role in building 
the skeleton of the G2 dynamical state space. Indeed, the earlier work of Hewitt and Wainwright [|29J 
provides some support for this expectation. 

In studying G2 cosmologies we expect to make use of insights into cosmological dynamics obtained 
from analysing SH cosmologies, for the following reasons. Gi cosmologies can be regarded as spatially 
inhomogeneous generalisations of SH models of all Bianchi isometry group types except Type- VIII and 
Type-IX, since, apart from these two cases, the G3 admits an Abelian Gi as a subgroup. In the language 
of dynamical systems the dynamical state space of SH cosmologies with an Abelian Gi subgroup is an 
invariant submanifold of the dynamical state space of Gi cosmologies. It thus follows that orbits in the 
Gi dynamical state space that are close to the SH submanifold will shadow orbits in that submanifold, 
thereby providing a link between Gi dynamics and SH dynamics. A further link is provided by the 
famous conjecture of BKL to the effect that near a cosmological initial singularity the EFE effectively 
reduce to ODE, i.e., the spatial derivatives have a negligible effect on the dynamics (see Bclinskii et al 
M, p656). In this asymptotic regime of near-Planckian order spacetime curvature it is plausible that SH 
dynamics will approximate Gi dynamics locally, i.e., along individual timelines. We shall refer to this 
conjecture as BKL II. 

As indicated above, we expect that SH dynamics will play a considerable role in determining the 
dynamics of Gi cosmologies, and that analogies with the SH case will be helpful. In two respects, 
however, the Gi problem differs considerably from the SH problem. Firstly, at any instant of time, the 
state of a Gi cosmology is described by a finite-dimensional dynamical state vector of functions of the 
spatial coordinate x. In other words, the dynamical state space of Gi cosmologies is a function space 
and, hence, is infinite- dimensional. The evolution of a Gi cosmology is thus described by an orbit in 

2 We refer to Hewitt and Wainwright p£|| for a dynamical formulation of perfect fluid G2 cosmologies using Hubble- 
normalised dependent variables. 

3 For details on the theory of FOSH evolution systems see, e.g., Courant and Hilbert [113] or Friedrich and Rendall |22| . 
4 See, e.g., van Elst et al U&. 
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this infinite-dimensional dynamical state space. The second difference lies in the complexity of the gauge 
problem, which we now describe. For SH cosmologies the G3 isometry group determines a geometrically 
preferred timelike congruence, namely the normal congruence to the G3-orbits, and hence there is a 
natural choice for the timelike vector field of the orthonormal frame. The remaining freedom in the 
choice of the orthonormal frame is a time-dependent rotation of the spatial frame vector fields, which we 
refer to as the gauge freedom. On the other hand, in a G2 cosmology there is a preferred timelike 2- space at 
each point that is orthogonal to the G2-orbits. Thus there is an infinite family of geometrically preferred 
timelike congruences and the gauge freedom in the choice of the orthonormal frame is correspondingly 
more complicated. There is also gauge freedom associated with the choice of the local coordinates. One 
of the goals of the present paper is to systematically discuss various gauge fixing options that arise for 
G2 cosmologies. 

The plan of the paper is as follows. In section |2| we derive the equation system that arises from the 
EFE and the matter equations. Working in the orthonormal frame formalism, we adapt the orthonormal 
frame to the G2-orbits and then simply specialise the general orthonormal frame relations in Ref. ||l9f to 
get the desired equation system in dimensional form. We then introduce the scale-invariant dependent 
variables and transform the equation system to dimensionless form. In section |3j we address the gauge 
problem and introduce four specific gauge choices, showing that our approach has the flexibility to 
incorporate all previous work. In section we discuss some features of the infinite-dimensional dynamical 
state space. In section || we give a simple geometrical representation of the past attractor as an invariant 
submanifold on the boundary of the infinite-dimensional dynamical state space. The nature of the past 
attractor illustrates the conjecture BKL II concerning cosmological initial singularities, and suggests the 
name " asymptotic silence''' to describe the dynamical behaviour of the gravitational field as one follows 
a family of timelines into the past towards the singularity. We conclude in section with a discussion 
of future research directions. Useful mathematical relations such as expressions for the scale- invariant 
components of the Weyl curvature tensor for G2 cosmologies and the propagation laws for the constraint 
equations have been gathered in an appendix. 

2 Framework and dynamical equation systems 

2.1 Dimensional equation system 

A cosmological model is a triple (Ai, g, u), where M. is a 4-dimensional manifold, g a Lorentzian 4-metric 
of signature ( — V + +), and u is the matter 4- velocity field. We will assume that the EFE are satisfied 
with the matter content being a perfect fluid with a linear barotropic equation of state, 

p(/5) = (7-l)£, 1<7<2. (1) 

The most important cases physically are radiation (7 = |) and dust (7 = 1). We will also include a 
non-zero cosmological constant A. Throughout our work we will employ units characterised by c = 1 = 
8ttG/c 2 . 

We assume that an Abelian G2 isometry group acts orthogonally transitively on spacelike 2-surfaccs 
(cf. Ref. J9J), and introduce a group- invariant orthonormal frame { e Q }, with e 2 and e 3 tangent to the G2- 
orbits. We regard the frame vector field eo as defining a future- directed timelike reference congruence. 
Since eo is orthogonal to the G2~orbits, it is hypersurface orthogonal, and hence is orthogonal to a 
locally defined family of spacelike 3-surfaces S:{t — const}. We introduce a set of symmetry-adapted 
local coordinates { t, x, y, z } that are tied to the frame vector fields e a in the sense that 

e = N- 1 d t , e 1 = e 1 1 d x , e A = e A B d x B, ,4,5 = 2,3, (2) 

where the coefficients are functions of the independent variables t and x only.H The only non-zero frame 
variables are thus given by 

N , e! 1 , e A B , (3) 

which yield the following non-zero connection variables: 

a, /3, ai, n + , er_, n x , er x , »-, «i, ^1 ; (4) 



5 In the terminology of Arnowitt, Deser and Misner Jfl, N is the lapse function, and we have chosen a zero shift vector 
field, N* = (i = 1, 2, 3). 
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their interrelation is given in the appendix. Here we have followed Ref. [[L9| in doing a (1 + 3)- 
decomposition of the connection variables. The variables a, /?, er_ and o~ x are related to the Hubble 
volume expansion rate H and the shear rate a a p of the timelike reference congruence eo according to 

a:=(H-2a+), p:=(H + a + ), (5) 

where <t+ is one of the components in the following decomposition of the symmetric tracefree shear rate 
tensor a a p: 



^4 



±(<T 22 +CT33) = -50U , V- ■= ^ (022 - 033) , &x ■■= ^<J 2 3 ■ (6) 

A consequence of this decomposition is that the shear rate scalar assumes the form a 1 := \[a a pa a ^ = 
3 {a\ + <r 2 + a x ) . We will use similar decompositions for the electric and magnetic Weyl curvature 
variables E a p and H a p, as given in the appendix. The "non-null-null" variables a and (3 (cf. Refs. 
p2| , na) turn out to be more convenient to use than H and ct+, since they are naturally adapted to the 
characteristic structure of the evolution equations that arise from the Ricci identities when the latter are 
applied to the timelike reference congruence eo (see Ref. p7|). The variables <X\, n + , n x and n_ describe 
the non-zero components of the purely spatial commutation functions a a and n a p, where 

n + := | (n 22 + n 33 ) , n_ := ^= (n 22 - n 33 ) , n x := -J= n 23 , (7) 

(see WE p3] for this type of decomposition of the spatial commutation functions) . Finally, the variable 
iii is the acceleration of the timelike reference congruence eo, while fii represents the rotational freedom 
of the spatial frame { e Q } in the (e 2 , e3)-plane. Setting £l\ to zero corresponds to the choice of a Fermi- 
propagated orthonormal frame {e a }. It should be pointed out that within the present framework the 
dependent variables 

{N, in, fix} (8) 

enter the evolution system as freely prescribable gauge source functions in the sense of Friedrich pl|, 
pl462 (see also section II.B of Ref. ||j). 

Since the G 2 isometry group acts orthogonally transitively, the 4- velocity vector field u of the perfect 
fluid is orthogonal to the G 2 -orbits, and hence has the form 

u = r(e + wei) , (9) 

where the Lorentz factor is T := (1 — w 2 ) -1 / 2 . It turns out to be useful to replace the matter energy 
density ft, in the fluid rest frame with |19| 

M=7T^tA, (10) 

(1 — v l ) 

where it is convenient to introduce the auxiliary quantities 

G±:=l±(7-1)« 2 . (11) 

Thus /i, which represents the matter energy density in the rest frame of eo, and v, which describes the 
matter fluid's peculiar velocity relative to the same frame, describe, for a given value of the equation of 
state parameter 7, the fluid degrees of freedom. 

The orthonormal frame version of the EFE and matter equations as given in Ref. |19| , when specialised 
to the orthogonally transitive Abelian G 2 case with the dependent variables presented above, takes the 
following form: 



Commutator equations 

Gauge fixing condition: 



0=(Cu)i:=N- 1 e l x d x N-in. (12) 



Evolution equations: 



N^dtd 1 = -ae/ (13) 

N^dte^ = -(f3 + V3a-)e 2 A -(V3a x +n 1 )e 3 A (14) 

N- x d t e 3 A = -(l3-V3o-^)e 3 A -(V3<T x -n 1 )e 2 A . (15) 
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Constraint equations: 

= (C com ) A 12 ~ {e 1 1 d x ~a 1 -V3n x )e 2 A -(n+-V3n_)e 3 A (16) 

= (C com ) A 3 i := (e 1 1 d x -a 1 +V3n x )e 3 A + (n+ + V3n-)e 2 A . (17) 

Einstein field equations and Jacobi identities 

Evolution equations: 

N- l d t a = -a 2 +/3 2 -3{a 2 _-n%+a 2 x -n 2 _)-al 

- \ 1 G^ l ^{l-v 2 ) + {e l 1 d x + Ui)ui (18) 

N~ l d t = -f /3 2 -§(tJ 2 +n 2 x + <t 2 +n 2_) _ i (2u 1 -ai)ai 

-ic;V[(7-i) + « 2 ] + 5 A ( 19 ) 

N- l d t ai = -/3(ui+a 1 )-3(n x cr_ -n_cr x )- §7G7_ 1 / u« (20) 

N~ dtn+ = — an + + 6 (cr_ n_ + (T x n x ) — (ei 1 9 X + ui) fii (21) 

N~ l d t a- + ei 1 9 K n x = - (a + 2/3) cr_ - 2 n+n_ - (ui - 2ai)n x - 2 fix ct x (22) 

AT - dtn x + ei d x a^ = — an x + 2cr x n+ — Ui cr_ + 2fiin_ (23) 

N~ l d t a x - e-i 1 d x n- = - (a + 2/3) cr x - 2n+n x + (wi - 2ai)n_ + 2Qi<t_ (24) 

A^~ 3tn_ — ei 9 X (7 X = — an_ + 2cr_ n.+ + ii\ a x — 2 Ox n x ■ (25) 

Constraint equations: 

= (G Gaufi ) := 2(2e 1 1 d x -3a 1 )a 1 -6(n 2 +n 2 _) + 2(2a + /3)/3-6(a 2 + a 2 ) 

- 2/i - 2A (26) 

= (Gcodacci)i := ei x d x (3 + ai (a-/3)- 3(n x <7_ -n_cr x ) - -jG^fxv . (27) 

Source Bianchi identities (Relativistic Euler equations) 

Evolution equations: 

— {N- 1 d t + -? r vei l d x )n + f 2 ex 1 d x v = - -J- /i [ a (1 + i; 2 ) + 2/3 + 2 (u x - oi) « ] (28) 

we/^Ju + Zaei 1 ^ = - -7^7-^ (1 - « 2 ) [ (2 - 7) <*v - 2 (7 - l)/3v 

+ G-u 1 + 2( 1 -l)a 1 v 2 ] , (29) 



/a 

/1 


{N- 1 


9t- 


Cx_|_Cx_ 


are 


(7- 


-1) 


(I-, 2 ) 2 , 



/2 : - ^2^ (1 - « 2 ) 2 . / 3 :=(37-4)-( 7 -l)(4-7)« 2 - (30) 

7Ur_ (jr_|_ 

Dynamical features exhibited by the dimensional evolution system for orthogonally transitive G 2 cos- 
mologies, that are independent of our later transformation to /3-normalised scale-invariant dependent 
variables, are the following: Eq. (|13[ ) evolves the only dynamically important frame variable (being part 
of the metric), Eqs. (|1§|) and ( |19| ) evolve the longitudinal components of the tensorial expansion rate 
of the timelike reference congruence eo, Eqs. (|2C| ) an d ( |2l| ) evolve a scalar and a non-tensorial spatial 
connection variable, respectively, Eqs. ( |22| ) - ( |25| ) provide the propagation laws for (transverse) gravita- 
tional waves, while, finally, the relativistic Euler equations ( |28| ) and ( p9| ) yield the propagation laws for 
(longitudinal) acoustic or pressure waves. Viewing the gauge source functions N, u\ and fii as arbitrarily 
prescribable real-valued functions of the independent variables t and x, the evolution system is already 
in FOSH format. The characteristic propagation velocities A relative to a family of observers comoving 
with the timelike reference congruence e aren 

A 1= 0, A 2 , 3 = ±l, A« = ^^,±fc^(7-l) 1/2 . (31) 



6 On characteristic propagation velocities, characteristic eigenfields, and related issues see, e.g., Refs. [1121] , 11 71 and |La|. 
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The right-propagating and left-propagating characteristic eigenfields associated with the non-zero A's are 
(for 1< 7 < 2)[] 

A 2 , 3 : 0_±n><) , (<r x T«-) , A 4 , 5 : j—, r(h 2 (-y,v)Thi{j,v)± _^ ) i/ 2 v ) > ( 32 ) 

where /ii(7, v) and ^2(7, v) are complicated expressions of their arguments that for v = have the limits 
/ii = and /i2 = 1, respectively. All of the A's are real-valued in the parameter range 1 < 7 < 2 of 
Eq. (|l|); this contains the dust case (7 = 1) and also the stiff fluid case (7 = 2). Note, however, that 
the former must be treated in terms of a modified version of the relativistic Euler equations, since for 
7=1 some of the coefficients in the principal part of the present version become zero. In summary, for 
our first-order dynamical formulation, the Cauchy initial value problem for the orthogonally transitive 
perfect fluid G2 cosmologies is well-posed in the range 1 < 7 < 2.^ 

The values of A 4i 5 reflect the anisotropic distortion of the sound characteristic 3-surfaces relative to the 
family of observers comoving with eo. This distortion may be viewed as a manifestation of the Doppler 
effect. In the limit v — ► 0, the magnitude | A 4j 5 | reduces to the isentropic speed of sound, c s = (7 — l) 1 / 2 . 
In the extreme cases 7=1 and 7 = 2 we obtain A4.5 = v and A4.5 = ± 1, respectively. 

Note that in this non-fluid-comoving form [even with the linear equation of state (0)] the principal part 
of the relativistic Euler equations is highly non-linear. This feature could lead to the formation of shocks 
in the fluid dynamical sector of the evolution system and should be kept in mind in a numerical analysis 
of the given equation system. The effective semi-linearity of the principal part of the gravitational field 
sector, on the other hand, is less likely to lead to numerical problems of this kind. For the latter jump 
discontinuities can be specified in the initial data for d x (cr- ± n x ) and d x (a x =F n_). 

The area density A of the G2-orbits plays a prominent role for G2 cosmologies. [] It is defined (up to 
a constant factor) by 

A 2 := (U a )(VbV b ) ~ (ZaV a ) 2 ■ (33) 

where £ and rj are two independent commuting spacelike Killing vector fields.Fj Expressed in terms of 
the coordinate components of the frame vector fields e^ tangent to the G2-orbits this becomes 

A' 1 = e 2 2 e 3 3 - e 2 3 e 3 2 . (34) 

The key equations for A, derivable from the commutator equations dlj) - (O) , are given by 

N-^ = 2(3, ei 1 ^ = -2a 1 , (35) 

i.e., j3 is the area expansion rate of the G2-orbits. The frame variables eA B , which play a subsidiary role 
as regards the dynamics of G2 cosmologies, are governed by Eqs. (|l4|) - (|l7|). Since these equations are 
decoupled from the remaining equations, we will not consider them further. 

2.2 Scale-invariant reduced equation system 

We will now introduce new dimensionless dependent variables that are invariant under arbitrary scale 
transformations. However, we will not normalise with the Hubble scalar H, as is usually done for SH 
models (here H is the volume expansion rate of the G3-orbits). Instead we will use the area expansion rate 
(3 of the G2~orbits, since this leads to significant mathematical simplifications for the resultant equation 
system in both its evolution and constraint part. We thus introduce /3-normalised frame, connection and 
curvature variables as follows: 



= (N-^e, 1 )/^ (36) 

= (tti, a x , a, <r_, n x , cr X) n_, n+, fii)//3 (37) 
= (^,A)/(3/3 2 ). (38) 



(AT 1 , iV) 
([/, A, (1-3E+), £_, N x , E x , 7V_, N + , R) 

(fi,fiA) 

7 Here we take the opportunity to correct for some sign errors in the expressions gisen in Refs. Jl 7| and |l8|. 

8 Clearly well-posedness is lost for any value of 7 in the range < 7 < 1 (cf . Ref . [B2| ) . 

9 The symbol we use to denote the area density, A, should neither be confused with the /3-normalised connection variable 
A to be introduced below, nor with the index "A" that takes the values 2 and 3. 

1 In terms of symmetry-adapted local coordinates x = y and x = z such that £ = d/dy and r] = d/dz, the area density 
is given by A = \/dct g^g, where 343, A, B = 2, 3, is the metric induced on the G2— orbits. 
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Note that we maintain the same notation that was introduced in WE |p3| for the iJ-normalised case. 
The two different normalisation procedures are linked through the relation 

#=(1 -£+)/?. (39) 

The above definition of E + in terms of a is motivated by the relation a = j3 — 3er+ , which follows from 
Eqs. ([|), and is equivalent to defining E + := er+//3. As a result of this definition, the /3-normalised shear 
rate scalar E 2 :— {a a pa a ^)/{<o(3 2 ) has the form E 2 = E+ + E 2 . + E 2 . Note that in the units we have 
chosen the matter variable v is already dimensionlcss. 



The dimensional equation system in subsection 2.1 leads to an equation system for the scale- invariant 



dependent variables ( p6| ) - (38). In order to make this change it is necessary to introduce the time and 
space rates of change of the normalisation factor [3. In analogy with /^-normalisation (see WE |5^]), we 
define variables q and r by 

M- l d t fi := -{q+l)/3 (40) 

= (C/j) := {E x l d x + r)[3 . (41) 

Here q plays the role of an "area deceleration parameter" , analogous to the usual "volume deceleration 
parameter", while r plays a role analogous to a "Hubble spatial gradient". Using Eqs. ( f40| ) and ( pi"] ) 
and the definitions (pq) - (pq), it is straightforward to transform the dimensional equation system to a 
/3-normalised dimensionless form. A key step is to use the evolution equation (fil3) for /3 and the Codacci 
constraint equation (|27| ) to express q and r, as defined above, in terms of the remaining scale-invariant 
dependent variables. The key result, which is essential for casting the scale-invariant evolution system 
into FOSH format, is that the expressions for q and r are purely algebraical . We refer to these equations 
as the defining equations for q and r [see Eqs. ( p4|) and ( |55| ) below]. The relation 

AT 1 d t r - E x l d x q = (q + 3E+) r - (r - U) (q + 1) + (q + 1) (C & ) (42) 

arises as an integrability condition for the decoupled /3-equations (Q) and (|4l|). 

Scale-invariant equation system 

Evolution system: 

Af-'dtE, 1 = ( g + 3E + )S 1 1 (43) 
3JV- 1 d t Y,+ = -3(g + 3E + )(l-E + ) + 6(E + + E 2 _ + E 2 ) 
+ | G; 1 [ (37 - 2) + (2 - 7 ) v 2 } n - 3 A 

- (E 1 1 d x -r + U -2A)U (44) 

M^dtA = (q + 3E+) A + (r - U) (45) 

Af- 1 8 t N + = (q + ST,+)N + + 6(T,_N_+E x N x )-(E 1 1 d x -r + U)R (46) 

A/^&Oa = 2{q+l)n A (47) 

AT' 1 d t E- + E-l 1 d x N x = (<? + 3E+-2)E_ -2A+A_ + (r-l> + 2,4)/y x -2i?E x (48) 

Af- l d t N x +£i 1 a E E_ = (<? + 3E + )/y x +2E X A+ + (r-L>)E_+2i?A_ (49) 

A/" _1 5 t E x -E x l d x N- = (g + 3E+-2)E x -2A+A^ X - (r - l> + 2A) A_ + 2i?E_ (50) 

A"" 1 d t N_ - E x l d x T, x = (g + 3E+)A_+2E_ N+-(r-U)T, x -2RN X . (51) 

^(AT^dt + ^vE^d^n + hE^d^ = 2-^/ 1 [^±(g+l)-i(l-3E + )(l + U 2 ) 
il G+ G+ 7 z 

-l + (r-U + A)v] (52) 

^0(A^- 1 9 t --4^z;S 1 1 9 x )i; + / 2 £; 1 1 a x n = 2 -A- Q (1 - v 2 ) [ { ^—^ (1 - « 2 ) r (53) 

7i G+G_ JiG_ 7 

- §(2- 7 )(l-3E+)v 

+ ( 7 -l)(l-Av)i;-|G_l7] , 



3 GAUGE CHOICES 9 

where /i, fa, f$ and G± are denned by Eqs. (|3(]) and (|ll|), respectively. 

Defining equations for q and r: 

q := ^ + ^(2U-A)A+l(j: 2 _+N 2 x +j:l+Nl) + l il '^ + V \ -ln A (54) 

r := -3AS + -3(JV x S--JV_E x )-|^Od. (55) 

Constraint equations: 

o = (Ccaufi) = fife - 1 + 2s + + s 2 . + s 2 + n + n A (56) 

= (C A ) = (£; 1 1 9 :E -2r)r! A , (57) 

where 

n h :=-%(Ei 1 d x -r)A + A* + Nl+Nl. (58) 

Gauge fixing condition: 

0=(Cjj) := AT' 1 E! 1 O x AT + (r - U) . (59) 



Supplementary equations 

The quantity (g+3S + ) occurs frequently in the scale-invariant equation system. Combining the definition 



(q + 3E+) = 2+{E 1 1 d x -r + il- 2A)A - § l - 7j (f - v 2 ) - 3H A . (60) 



of q given in Eq. (54) with the Gaufi contraint equation ( J56| ) solved for E+, one can express this quantity 

by 

(2_7) 
G 4 

In terms of our scale-invariant dependent variables, the area density A of the G2~orbits satisfies the 
relations 

A- 1 Af- 1 d t A = 2, A~ 1 E 1 1 d x A = -2A. (61) 

Combining the two, the magnitude of the spacetime gradient V ' a A is 

{V a A){V a A) = -Af3 2 (\-A 2 )A 2 , (62) 

so V a A is timelike for ^4 2 < 1. 

3 Gauge choices 

In this section we discuss the gauge problem. 
3.1 Gauge freedom 



The scale-invariant equation system in subsection 2J2 contains evolution equations for the dependent 
variables 

{£i\ S+, A, N +1 J2 Al 2_, JV X , E x , JV-, n,v} , (63) 

but not for the gauge source functions 

{M, U, R}, (64) 

and thus does not uniquely determine the evolution of the G2 cosmologies. The reason for this deficiency 
is that the orthonormal frame { e a } and the local coordinates { t, x } were not specified uniquely in 



subsection 2.1. We now summarise the remaining freedom, which we refer to as the gauge freedom. 

(i) Choice of timelike reference congruence eo. 

The gauge freedom is a position-dependent boost 

t)- T (l ?)(;). -=7^' -«*•'>• < 65 » 

in the timelike 2-spaces orthogonal to the G2~orbits. 
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(ii) Choice of local time and space coordinates t and x. 
The gauge freedom is the coordinate reparametrisation 

i=i(t) , x = x{x) . (66) 

(iii) Choice of spatial frame vector fields e 2 and e3 . 

The gauge freedom is a position-dependent rotation 

e 2 \ _ / cos ip shop \ ( e 2 
e 3 J \ -simp cos<^ J \ e 3 

in the spacelike 2-spaces tangent to the G 2 -orbits. 



tp = p>(t,x) , (67) 



We say that (i) and (ii), which refer to the freedom associated with the preferred timelike 2-space, 
constitute the temporal gauge freedom, and that (iii) represents the spatial gauge freedom. 

Table ti\ shows possible ways of fixing the spatial gauge by requiring one frame vector field or a 
combination of frame vector fields to be parallel to a Killing vector field. Each of these sets of conditions 
is preserved under evolution and under a boost. These choices are essentially all equivalent. We will 



Spatial gauge condition Frame vector field 

parallel to a KVF 

N+ - V3N- = = R + n/JKx &2 

N+ + V37V_ =0 = R- V3S X e 3 

N + + V3N X = = R + \/3£_ e 2 - e 3 

N+ - \Z3JVx =Q = R- V3E_ e 2 + e 3 

Table 1: Spatial gauge conditions for aligning a combination of the frame vector fields e^ with a KVF. 
routinely make the first choice, namely 

N+ = VSN-, R = -V3E X . (68) 



With this choice the evolution equation (Utj) becomes identical to Eq. ( |51| ) , and thus can be omitted from 
the full scale-invariant equation system. Other interesting choices for fixing the spatial gauge do exist, 
however, such as a Fermi-propagated frame, for which R = 0. 

3.2 Fixing the temporal gauge 

Within the present scale-invariant formulation of the dynamics of orthogonally transitive G 2 cosmologies 
with perfect fluid matter source we will fix the temporal gauge by adapting the evolution of the gauge 
source function U to the following four geometrical features, listed in order of subsequent discussion. 

(i) Adapt the evolution to a family of null characteristic 3-surfaces. 

(ii) Adapt the evolution to the integral curves determined by the spacetime gradient of the area density 
of the G 2 -orbits, V a A 

(iii) Adapt the evolution to the family of fluid sound characteristic 3-surfaces. 

(iv) Adapt the evolution to zero-velocity characteristic 3-surfaces associated with a family of freely- 
falling observers. 

The idea is to specialise eo in such a way that cither M^ 1 dtU or U itself is determined in terms of 
the other dependent variables. Then JV is determined from Eq. ( p9| ) up to an arbitrary dimensionless 
multiplicative function f(t). We then use a reparametrisation of t to choose 

f{t) = e ct , (69) 
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where C is an arbitrary constant. This coordinate choice leads to an autonomous differential equation 
for A/" which we include in the evolution system, giving a fully determined autonomous scale-invariant 
equation system. 

It should be pointed out that apart from the above four choices of temporal gauge other interesting 
possibilities such as, e.g., a constant area expansion rate gauge, where r = 0<!=)>/3 = Pit), do exist. 

3.2.1 Null cone gauge 

The first choice of gauge, which we call the null cone gauge, is motivated by the identity 

Af- 1 d t (r -U)- E, 1 d x (q + 3E+) = - AT 1 E x x d t d x [ln(AfE 1 1 )]+Af- 1 d^) , (70) 

which follows from combining Eqs. (B9[) and (f43|). It suggests that we impose the condition 

Q = M- 1 d t {r-U)-E l 1 d x (q + -i^+) (71) 

on A/" -1 dtU. It follows immediately from these two relations and Eq. (E3) that 

N = ^M (72) 

El 

and 

\ AT 1 d t U + E^ d x Z+ = \{q + 3E+) U - | (r - U) (1 - 3£+) , (73) 



provided that the gauge fixing condition ( |59|) propagates along eo according to Eq. (188) in the appendix. 
We now use the i-reparametrisation ( p6[ ) to set /(£) = e Cnct , where C nc is a constant. Equation ( |72; ) then 
gives 

MEi l =e a ^*g{x) . (74) 

On differentiating Eq. ( |74| ) and using Eq. (Q) , we obtain an evolution equation for Af that reads 

Af- 1 d t Af=-(q + 3J: + )Af + C nc . (75) 

Note that in the null cone gauge Eqs. ( p4[ ) and ( J73J ) form the (£+, t7)-branch of an autonomous evolution 
system in FOSH format. The associated characteristic propagation velocities are A = ± 1. 

Choosing the null cone gauge permits one to introduce the familiar conformal coordinates { t, x } in 
the timelike 2-spaces orthogonal to the G2~orbits, although we do not find it convenient to make this 
choice in general. Referring to Eq. (172), one can use the coordinate reparametrisation ( pq ) to set f(t) = 1 
and g(x) = 1, so that AfE] 1 = 1. It follows from Eq. (Bq) that Nei 1 = 1, which implies, using Eqs. (0), 
that the line element in the timelike 2-spaces orthogonal to the G2~orbits has the form 

^ds 2 = N 2 (- dt 2 + dx 2 ) . (76) 

Conformal coordinates have been frequently used in the analytical study of vacuum Gi cosmologies, and 
in the derivation of exact solutions, both for vacuum and for perfect fluid models. Selected references 
from the literature are Gowdy J24J, g5[ , Liang |35|] , Isenberg and Moncrief J33[ , Hubner |32J , Kichenassamy 
and Rendall [|34| , Senovilla and Vera J47[ and Anguige 0] . 

3.2.2 Area gauges 

The separable area gauge is determined by imposing the condition 

0=(r-U), (77) 



which determines U algebraically through Eq. (p5|). There is thus no need to determine an evolution 
equation for U. It follows immediately from the gauge fixing condition (|5|) that Af = f(t). We now use 
the i-reparametrisation (|6q ) to set /(£) = A/"o, a constant, i.e., 

N = N . (78) 
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In this case the evolution equation for Af is trivial, i.e., 

A/"" 1 d t Af = . (79) 

It follows from Eqs. (pi]) and (p8|) that the area density has the form 

A = e 2 Q e 2Mot m(x) , (80) 

where here and throughout £q denotes the unit of the physical dimension [length], a nd m (x) is a positive 
function of x. The gauge fixing condition (|77j) propagates along eo according to Eq. ( |l89| ) in the appendix 
subject to an auxiliary equation for A/g" dtJJ. Note that the separable area gauge does not in general 
yield an evolution system in FOSH format. 

For the class of Gi cosmologies in which the spacetime gradient V a A is timelike, we can strengthen 
the separable area gauge condition (I77J) by requiring in addition that 



A = 0, (81) 

which we achieve by choosing eo to be parallel to V a A. It follows from Eq. ( pl| ) that d x A — 0, and Eq. 
( p0| ) reduces to 

A = £ 2 e 2Not . (82) 



This defines the so-called area time coordinate. Observe that condition ( jSjJ ) is invariant, by virtue of 
Eqs. ( (45J ) and (77). We shall refer to the gauge choices ([77]) and (81) as the timelike area gauge. We note 



that in this case, with S + and U algebraically determined in terms of the other scale-invariant dependent 
variables from, respectively, the Gaufi constraint equation ( pq ) and Eqs. ([77]) and (pq), the evolution 
system becomes unconstrained when £1a = and does assume FOSH format. We give the resultant 



equation system in subsection 1.4. Selected references from the literature using the timelike area gauge 



are Berger and Moncrief M, Hern and Stewart M and Rendall and Weaver |43j]. 

3.2.3 Fluid-comoving gauge 

The fluid-comoving gauge is determined by choosing eo to be equal to the fluid 4-velocity field u. By 
virtue of Eq. (|9|), this choice is equivalent to imposing the condition 

= v . (83) 

The evolution equations (|5^) and ([33]) for tf, and v now reduce to 

N- l d t n = 2[( ff + 3E+-2) + |(2-7)(l-E + )]n (84) 

= (C) := [(7-l)(E 1 1 d x -2r)+ 7 U\ fi. (85) 

The evolution equation for U (which is now identified with the /3-normalised fluid acceleration) results 
from demanding that the new constraint equation ( p5j ) propagates along u. This leads to 

i AT 1 d t U + ( 7 - 1) E x x d x V+ = \{q + 3S+) U - ( 7 - 1) (1 - £+) {r - U) . (86) 



It then follows that the gauge fixing condition ( |59| ) propagates along u according to Eq. (190) in the 
appendix. Furthermore, it follows from Eqs. (JS5|) and (p9) thatFH 



»2o2 r ,\-h- 1 )/^f 



Af = f{t)(^){^ 2 Vt) n 1,n . (87) 

We now use the t-reparametrisation ( pq ) to set fit) = e Cfc *. On differentiating Eq. (pTT) and using Eqs. 
( p0| ) and (p3) , we obtain an evolution equation for Af that reads 

Af- 1 d t Af=-[(q + 3S+ - 2) + 3 (2 - 7) (1 - £+) ] Af + C tc ■ 



11 This equation is a scale-invariant form of the well-known dimensional relation N = a(t) [i ' 7 1 ^' 7 for perfect fluid 
models with equation of state pip) = (7 — 1) /i in fluid-comoving gauge; see, e.g., Ref. [BfJ. 
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In order to obtain an evolution system in FOSH format, we need to multiply Eq. ( |44| ) by a factor of 
(7 — 1) and write it in the form 

( 1 -\){'iN- 1 dj: + +E 1 1 d x U) = - (7 - 1) [ (q + 3£+)(l-3£+) + 2g- 6 (£?. + £*) 

- (r - U + 2A) U - § (3 7 - 2) Q + 3f2 A ] ■ (89) 



When we adjoin Eqs. (]8q), fl88| ) and (p9J) to the full scale- invariant equation system, with the fluid 
evolution equations (p2l) and (p3\) replaced by Eq. (p4|) and the new constraint equation (J85l) , and Eq. (E3) 



replaced by Eq. (89), we obtain a new evolution system that has FOSH format, with the fluid dynamical 
sector being shifted from the (O, i>)-branch to the (£+, E/)-branch, with characteristic velocities given by 
A 4 , 5 = ± (7 — l) 1 / 2 . The family of sound characteristic 3-surfaces thus becomes symmetrically embedded 
inside the family of null characteristic 3-surfaces. In the case of dust, 7 = 1 <=> U = 0, Eq. (g9|) does 
apply without the common factor (7 — 1). 

When doing numerical experiments in the present framework, the fluid-comoving gauge will have the 
advantage that, in view of the linear equation of state (|lj), the effective semi-linearity of the principal 
part of Eqs. (g9[) and (pq ) will prevent the development of shocks in the fluid dynamical sector of the 
evolution system. Hence, only the propagation of so-called contact discontinuities is possible in both the 
gravitational field and the fluid dynamical sectors. Another advantage of this gauge is that it makes 
direct physical interpretation possible in terms of kinematical fluid quantities. 

Examples of references employing the fluid-comoving gauge are Eardley et al [H, Wainwright and 



Goode |54] and Ruiz and Senovilla 45 



3.2.4 Synchronous gauge 

The synchronous gauge is determined by choosing eo to be a timelike reference congruence that is geodesic , 
i.e., we impose the condition 

= U . (90) 



It follows from the gauge fixing condition (59) and Eq. (|4l|) that 



N = fit) (4/3) • (91) 



We now use the i-reparametrisation ( pq ) to set f(t) = e Csync *. On differentiating Eq. (J9l]) and using Eq. 
( ffO| ) , we obtain an evolution equation for M that reads 

N- l d t N=-{q + l)N + C sync . (92) 



The gauge fixing condition fl59| ) presently propagates along e according to Eq. ( |191| ) in the appendix. 
When we adjoin Eq. (p2n to the full scale-invariant evolution system, simplified using Eq. (pfl), we again 
obtain FOSH format. In a more general context, the synchronous gauge has been made prominent in 
particular by the work of BKL J5L pf . 

4 Scale-invariant dynamical state space 

The description of the dynamics of Gi cosmologies is complicated by the fact that the scale-invariant 
dynamical state vector, and hence the structure of the scale-invariant dynamical state space, depends 
on the choice of gauge. In this section we discuss this issue, and we explain which properties of the 
scale- invariant equation system and of the dynamical state space are independent of the choice of gauge. 
For simplicity, in the present discussion we set the cosmological constant to zero, J7a = 0. 

4.1 Overview 

We assume that the spatial gauge has been fixed according to Eq. (pq). Once we choose a specific 



temporal gauge, the equation system derived in subsection 2/1 gives an explicit set of evolution and 
constraint equations for a finite-dimensional dynamical state vector X. These equations can be written 
concisely in the following form, where the FOSH nature of the evolution part is indicated by the fact that 
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the coefficient matrices A(X) and B(X) are symmetric, with A(X) being positive definite. 

Evolution system: 

A(X) a X + B(X) diX = F(X) , 



(93) 



and 

Constraint equations: 

0=C(X,fl]X), (94) 

where do := A/" -1 <9t and 3i := Ei 1 d x . The dynamical state vector X depends on the choice of temporal 
gauge as follows: 

X = X g © X w , (95) 

where X g is the temporal gauge-dependent part, and 

X W = (£_,7V X ,£ X , N^f , (96) 

the temporal gauge-independent part describing the dynamical degrees of freedom in the gravitational 
field. The latter amount to four arbitrary real-valued functions of x one can specify initially. 



Temporal 
gauge 


x g 


Constraint 
equations 


Initially freely 
specifiable 
functions 


No. of gauge 

degrees 

of freedom 


Fluid-comoving 


(Ar,E 1 1 ,z + ,u ) A,n) T 


(C[>), (Ccaufj), (C v ) 


{Ei\A,£l) T 


1 


Timelike area 


(Ei l ,n,vy 


none 


(E 1 \n, v y 


1 


Separable area 


(SiSE+^fi,*) 7 ' 


(Ccaufi) 


(E x \A,n,vy 


2 


Synchronous 


(A^AE-nA^r 


(Cjj), (C G aufi) 


{E-l^A^v) 1 


2 


Null cone 


(AA,£i\ £+,£/, A ^,f) T 


(Cjj), (Ccaufi) 


(£i 1 ,iM,n,v) T 


3 



Table 2: Form of X g and number of gauge degrees of freedom for different temporal gauge choices. 

In Tab. [2| we show the number of gauge degrees of freedom that remain after choosing a specific temporal 
gauge. This number is arrived at as follows: 

No. (initially freely specifiable functions) = Dim(X g ) — No. (constraint equations) 

No. (gauge degrees of freedom) = No. (initially freely specifiable functions) — 2 . 

We now discuss how the remaining gauge degrees of freedom arise. Once a temporal gauge has been 
chosen, each set of initial conditions 

cX = X(t ,x) 

that satisfies the constraint equations determines a unique solution of the evolution equations. Because 
of the remaining gauge freedom, different initial conditions do not necessarily lead to physically distinct 
solutions. We can use the remaining gauge freedom to simplify the initial conditions as follows. Firstly, 
the £-reparametrisation d6j1) can be used to set 



.Ei 1 (to, a;) = 1 for all x 

in each gauge, thereby eliminating one gauge degree of freedom. At this stage there is no remaining gauge 
freedom in the fluid-comoving and timelike area gauges. The remaining gauge freedom in the three other 
temporal gauges is a boost (65) with velocity w — w(t,x) that preserves the appropriate gauge fixing 
condition. The requirement that the gauge fixing condition be preserved leads to a propagation equation 
for w(t,x), of first order in time for the separable area and synchronous gauges and of second order in 
time for the null cone gauge. Thus, for the first two gauges one has the freedom to perform a boost at 
the initial time with w(to,x) arbitrary, while in the null cone gauge both w(to,x) and dtw(to,x) can be 
chosen arbitrarily. One can use this restricted boost with w(to,x) arbitrary to set v(to,x) — 0, and in 
the case of the null cone gauge one can use the arbitrary function dtw(to,x) to also set U(to,x) = 0. 
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In other words, in the separable area and synchronous gauges one can, without loss of generality, use 
the initial condition v(to, %) = 0, and in the null cone gauge one can likewise, without loss of generality, 
use the initial conditions v(to,x) = and U(to,x) — 0. Thus, in each gauge, X g is specified initially by 



giving two arbitrary real- valued functions of x, which, with X w in Eq. (96), gives a total of six dynamical 
degrees of freedom. 

4.2 Familiar solutions as invariant submanifolds 

The Gi cosmologies contain a rich variety of familiar classes of solutions as special cases. In this subsection 
we indicate how these classes of solutions arise as invariant submanifolds in the scale-invariant dynamical 
state space. In each case the related equation system can be obtained by specialising the general scale- 
invariant equation system and making an appropriate choice of gauge. 

4.2.1 Vacuum Gi cosmologies 

These solutions are described by the subset 

= O, (97) 

which is invariant since O = implies 3(0 = by Eq. (p2[). If the spacetime gradient V a A is timelike, one 
can use the timelike area gauge. The resulting unconstrained evolution system, which has FOSH format, 



is given in subsection 4.4 [Eqs. (123) - (125)]. If the spatial topology is T 3 , these equations describe the 
Gowdy vacuum spacetimes that can contain gravitational radiation with two polarisation states J24, E5| . 
It should be noted that the evolution equation (|3|) for v is singular on the vacuum boundary O = 
due to the fact that if one solves for dtV, one obtains the singular term <9 x O/0. This fact means that care 
has to be taken in taking limits as O — > 0, unless one is working in the fluid-comoving gauge, in which 
case this problem does not arise. 

4.2.2 Diagonal G2 cosmologies 

The orthogonally transitive G2 cosmologies have in general four dynamical degrees of freedom in the 
gravitational field, of two different polarisation states, that are associated with the null characteristic 
eigenfields (£_ ± iV x ) and (£ x T N-). With the spatial gauge choice (p8[), it follows that the conditions 



= S x = 7V_ (98) 

define an invariant submanifold which corresponds to Gi cosmologies with one possible polarisation state 
only. We shall refer to this class of solutions as diagonal G2 cosmologies, because for them the line 
element can be written in diagonal form (since both Killing vector fields are hypersurface orthogonal; cf. 
WE §§). 

4.2.3 Plane symmetrical Gi cosmologies 

Specialising further, we can eliminate both polarisation states by considering the invariant submanifold 

= S X = N_ =£_ =iV x , (99) 

which describes the class of plane symmetrical Gi cosmologies (the isometry group here is a G3 acting 
multiply-transitively on flat spacelike 2-surfaces). These solutions are the plane symmetrical analogues of 
the well-known spherically symmetrical Lemaitre-Tolman-Bondi models, in general with non-zero fluid 
pressure (see Stewart and Ellis Q an d Eardley et al Q ) . When 7 = 1 the evolution system reduces to 
a set of ODE. 

4.2.4 Self-similar G2 cosmologies 

It is of interest to consider the Gi cosmologies that correspond to the equilibrium points (i.e., fixed points) 
of the evolution system (B3J), that are defined by the condition 

= d X. (100) 
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This condition means that the dynamical state vector X is constant on those timelike 3-surfaces whose 
spacelike normal congruence is ei. It follows that these 3-surfaces are the orbits of a 3-parameter homo- 
thety group H3, i.e., solutions of this kind are self-similar. It is important to note that the condition ( |100[ ) 
should be imposed before specifying the temporal gauge, since it uniquely fixes the gauge by specifying 
the timelike frame vector field eo. Indeed, Eq. ( |100| ) implies that both the separable area gauge condition 
( f77| ) and the null cone gauge condition (^Tj) are satisfied. 

Under these conditions, the evolution system (J93|) reduces to a set of ODE that govern the spatial 



dependence of the models. In other words, the condition (100) defines a finite-dimensional submanifold of 
the infinite-dimensional dynamical state space. Specialising further, the condition v = defines a smaller 
invariant set of solutions consisting of self-similar models whose fluid 4-velocity field u is tangent to the 
i?3-orbits. These solutions, which we shall refer to as fluid- aligned self-similar Gi cosmologies, have been 
analysed qualitatively in some detail by Hewitt et al |2£j, |3l], [27J . They are of interest as potential future 
asymptotic states for more general G2 cosmologies. 

4.2.5 Spatially homogeneous Gi cosmologies 

The conditions 

= diX, = U , (101) 

define a finite-dimensional invariant submanifold of the infinite-dimensional dynamical state space cor- 
responding to Gi SH models, which admit a G3 isometry group acting transitively on the spacelike 



3-surfaces orthogonal to eo.fjj As with Eq. (|100j ), the condition (101) should be imposed before speci- 



fying the tempor al g auge, since it likewise fixes the gauge by specifying the timelike frame vector field 



eo. Indeed, Eq. (101) implies that both the separable area gauge condition ( |77| ) and the null cone gauge 
condition (^Tj) are satisfied. 

Under these conditions, the evolution system ( |93] ) reduces to a set of ODE which determines the 
dynamical evolution of the models. Since d\K = 0, the _Ei 1 -equation (J43|) decouples from the full 
system. Thus, if one is only interested in the evolution of Gi SH models, all relevant information is 
given by the remaining equations, which are analogous to the equation systems studied in WE, but with 
iJ-normalisation replaced by /3-normalisation. However, since we are interested in how the G2 SH models 
are related to the Gi cosmologies, it is necessary to retain the i?i 1 -equation. Specialising further, the 
condition v — defines a smaller invariant set of solutions, the so-called non-tilted SH models, in which 
the fluid 4-velocity field u is orthogonal to the G3-orbits. 

The equilibrium points of SH dynamics, i.e., cosmologies that admit an H4 acting transitively on 
spacetime, play an important role in the Gi dynamical state space. There are two main subclasses. 
Firstly, those equilibrium points that satisfy (g + 3S + ) ^ must satisfy Ei 1 — 0, on account of Eq. (fT~ 



They are thus constrained to lie in the unphysical boundary E\- = (see subsection 4.3), and hence 



can potentially affect the Gi dynamics near the cosmological initial singularity. The most important 



examples are the Kasner equilibrium set (see subsections 4.3 and 5.2) and the flat FL equilibrium point 



(see subsection 5.4). Secondly, those SH equilibrium points that satisfy (q + 3£+) = lie in the physical 
part of the dynamical state space and hence can potentially affect the evolution at late times. The most 
important of these are the so-called plane-wave equilibrium points (see WE [631, Ch. 9). 

In the present formulation it is possible to solve globally for E + from the Gaufi constraint equation 
(p6|). So, for example, in the vacuum case one thus automatically obtains a reduced dynamical system 
whose dimension is equal to the number of dynamical degrees of freedomJ 13 | 

It is worth noting that if one introduces the standard Fermi-propagated diagonal frame for SH models 
of class A, which is not the default frame choice in our formulation, the Kasner set is represented by a 
parabola given by 2S + + £_ = 1, while the Type-II vacuum solutions are straight lines. 

4.3 Unphysical boundary 

The evolution equation (f43| ) for the scale-invariant frame variable E\ 1 shows that the set E\ 1 = defines 
an invariant submanifold in an arbitrary gauge. This invariant submanifold divides the dynamical state 

12 All perfect fluid SH cosmologies of Bianchi Type-I to Type-VIIjj, apart from the exceptional Type-VI_i/g, admit an 
Abelian G2 subgroup which acts orthogonally transitively, and are hence included. 

13 Thus there exists no drawback with the orthonormal frame approach as recently indicated by Szydlowski and Demarct 
|51| (of course, a /3-normalisation can also be done in the Type- VIII and Type-IX cases). 
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space into two disjoint invariant submanifolds given by Ei 1 > and Ei 1 < 0. The full scale-invariant 
equation system is, however, invariant under the discrete symmetry 

(x, Si 1 ) — » {-x, -E^). (102) 

The two invariant submanifolds are thus physically equivalent, and without loss of generality we can 
restrict our considerations to the case 

£a x >0. (103) 

The set E\ x — corresponds to unphysical states for which the area expansion rate (3 diverges (/?—►+ oo), 
typically leading to a spacetime singularity. We shall refer to the invariant submanifold E\ l = as the 
unphysical boundary of the infinite-dimensional dynamical state space. It is significant that the evolution 
system is well-defined on the unphysical boundary E\ l =0. Indeed Eq. (|93| ) reduces to 

A(X)d X = F(X) , (104) 



i.e., a system of ODE, on the unphysical boundary. It is important to note that the solutions of Eq. (IO4), 
regarded as solutions of the full evolution system (p3|), have arbitrary x- dependence. ,Fj These solutions 
thus do not, in general, correspond to solutions of the EFE, and in this sense they are unphysical. 
Nevertheless, they do play a significant role in the evolution of G2 cosmologies. The key point is that 
if an orbit in the physical part of the dynamical state space with E\ l > approaches the unphysical 
boundary as t — > — 00, then it will shadow orbits in this boundary, i.e., the dynamics in the unphysical 
boundary will determine the asymptotic dynamics of Gi cosmologies that are solutions to the EFE and 
are thus regarded physical. In the unphysical boundary there is a hierarchy of invariant submanifolds that 
influence the asymptotic dynamics of Gi cosmologies. Firstly note that on the unphysical boundary the 
gauge fixing condition (|59|) reduces to (r — U) — 0, which, in conjunction with the integrability condition 
([44), implies that 

Af- 1 d t U= (q + 3Y, + )U . (105) 

It follows that, in an arbitrary gauge, the condition [7 = defines an invariant submanifold in the 
unphysical boundary. The evolution equations for the invariant submanifold = E\ 1 = U are precisely 



the evolution equations for SH models, as follows from subsection 4.2. 

Within the invariant submanifold = E\ l = U, the vacuum subset Q = is invariant, and within 
this set is the Kasner invariant set, defined by 

= A = N x = N_ . (106) 



These conditions imply, on account of the Gaufi constraint equation (56) and Eqs. ( |58| ) and (pOJ), that 

( 9 + 3E+)=2, (107) 

and 

(108) 



2E+ + Ef. 


+ K = 1 ■ 


The remaining evolution equations are 




M- X d t ^ + = 





Af- 1 <9 f £_ = 


2V3£^ 


AT-^Ex = 


-2\/3S_S 



(109) 
(110) 
(111) 

where the dependence of S + , £_ and S x on the local coordinate x is unrestricted. The Kasner equilibrium 
points are given by 

E x =0, 2S+ + El = l > (112) 

where £_ is an arbitrary function of x. These conditions define a Kasner parabola /C. Intuitively speaking, 
the orbits in the Kasner invariant set, including the Kasner equilibrium points, describe a G2 cosmology 
with the evolution of a Kasner vacuum solution of the EFE, but with unrestricted ^-dependence. 

14 An important example of such a solution is a Kasner metric whose Kasner exponents, instead of being constants, depend 
on the local coordinate x. 
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4.4 Timelike area gauge 

In this subsection we give the evolution system in the timelike area gauge, which was introduced in 
subsection 3.2.2J . For simplicity, we assume that the cosmological constant is zero, ^a = 0. This gauge 



has the advantage of leading to an unconstrained evolution system in FOSH format, as follows. 

No 1 d t E 1 l = (q + S&jEh 1 (113) 

7VV 1 <9 t E_ + E x l d x N x = (g + 3E+-2)E_+2\/3E^ - 2V^iV£ (114) 

TVq- 1 d t N x + V X S_ = (<? + 3£ + )7V x (115) 

N^ 1 dt^x - Ei 1 d x N- = (g + 3S+-2-2\/3S_)S x -2y/3N x N- (116) 

No 1 d t N- - Ex 1 d x T, x = (g + 3£+ + 2V3£_)JV_+2\/3£ x iV x (117) 

^{M^dt + ^vE^d^n + hE^d^ = 2-^/i (118) 

\l Cjt_}_ Ct_j_ 

x [^± (g+ l)-i(l-3E + )(l + W 2 )-l] 

7 

^O(AA - 1 a t --^«^i 1 ^)« + / 2 i?i 1 a^ = -J^£l{l-v 2 )[ { ^—^G + U 
h G+G- jiG- 7 

+ (2- 7 )(l-3E + )»-2(7-l)«] . (119) 



where 



3(2-7),, .. 2 x 



( 9 + 3E+) = 2- |^—^(l-^)0. (120) 

G+ 

The auxiliary variables E + and t/ are obtained from the Gaufi constraint equation (Bq) and Eqs. ( [77] ) 
and ( Ba ) as 

E+ = \ (1 - E 2 . - N\ - E 2 - Nl - Q) (121) 

U = r = -3(JV x E_-7V_E x )-f-^ftw . (122) 

G+ 



Note that in the present case we have from Eqs. fll20| ) and (121) that q > i which, on account of Eq. 
([40|), guarantees that (3 is a monotone function. 

AAA Gowdy vacuum spacetimes 



The vacuum subcase of Eqs. (|113j ) - ( 120 ) describes amongst others the Gowdy spacetimes with spatial 



topology T 3 |2J, |25| . In characteristic normal form these equations can be written asQ 

JV^dtEi 1 = 2E X X (123) 

{Afo 1 d t ±E 1 1 d x )(Z_±N x ) = (E_±7V x )-(E_TA r x) + 2x/3(E xT iV_)(E x ±iV_) (124) 

(Afo 1 d t ±E 1 1 d ;€ )(Z x TN_) = (E xT A r -)-[l + 2\/3(E_±7V x )](E x ±iV_) . (125) 



Note that in the present case we can use a reparametrisation ( p6|) of a; to set d x E\ Y — 0. We use this 
representation to exemplify the following three aspects, which hold for the whole class of Gi cosmologies 
and, in suitably generalised form, indeed for any general cosmological model. 



(i) As the source terms on the RHS of the gravitational field equations (124) and ( |125| ) (and E\ l ) must 
be continuous for the PDE system to be well-defined in the ordinary sense, so are the four charac- 
teristic first derivatives on the LHS. This implies that there are four unrestricted first derivatives 
given by (A^ -1 d t =p Ex 1 d x ) (E_ ± N x ) and (A/" _1 d t T E x x d x ) (E x =p N_), which can be thus in- 
terpreted as the arbitrary information (four free real-valued functions) that gravitational radiation 
can propagate. 



15 On the characteristic normal form of a FOSH evolution system, see Ref. |12|. 
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(ii) Spacetimes with = (£ x — 7V_) = (S x + JV_) do form an invariant submanifold of the dynamical 
state space; as mentioned before they correspond to the diagonal subcase for which the dynamical 
degrees of freedom in the gravitational field are "polarised" . 

(iii) Right-propagating and left-propagating characteristic eigenfields of the gravitational field do not 
form invariant submanifolds of the dynamical state space; they cannot in general be separated from 
each other. Referring to the Weyl curvature components listed in the appendix, this reflects the 
fact that a typical G2 cosmology (and any general cosmological model) is of algebraic Petrov type I. 

5 Past asymptotics and the past attractor 

In this section we discuss the asymptotic evolution of the class of orthogonally transitive Gi cosmologies 
near the cosmological initial singularity. We present evidence to support the claim that the past attractor 
lies on the unphysical boundary, and is a subset of the Kasner parabola JC. This analysis leads to a 
discussion of the notion of asymptotic silence, and the related conjecture BKL II. Because of orthogonal 
transitivity of the G2 isometry group, we do not expect to find a past attractor of oscillatory nature. We 
also linearise the evolution equations about the flat FL equilibrium point, and relate the results to the 
concept of an isotropic initial singularity. 



5.1 Past attraction to the unphysical boundary 



We use the tim elike area gauge characterised by E qs. (|77|), J78[ ) and (plh, and the as socia ted evolution 
system (113) - (12C). The evolution equation ( 113 ) for Ei 1 , in conjunction with Eq. (120), reads 



Mo 1 d t E l l = 



2-l^d 



v 2 )n 



Et 1 . 



(126) 



For vacuum models, i.e., il = 0, we can solve this ODE, obtaining Ei 1 — b(x) exp(2JVot), which implies 

(127) 



lim Ex 1 = 

t — >— 00 



This equation will also hold for perfect fluid models which satisfy the requirement^] 



lim O = 

t — > — OO 



(128) 



If the orbit of a G2 cosmology satisfies Eq. (127), the orbit will approach the unphysical boundary, 
and then we expect that it will shadow orbits in the boundary, which are described by the system of 
ODE obtained from Eq. @ by setting E x 



0. On the unphysical boundary U satisfies the evolution 
equation (105), which is of the same form as Eq. (126). It follows that along orbits in the unphysical 
boundary limt^-oo [/ = 0, and so we expect that typical orbits will approach the invariant submanifold 
U. As mentioned in subsection 4.3, the evolution system in this invariant submanifold is 



precisely the evolution system for SH models. We thus expect that SH dynamics will approximate Gi 
dynamics asymptotically as t — ► — 00. These heuristic considerations suggest that we should consider the 
Kasner parabola K in order to localise a possible past attractor. 



5.2 Linearisation about the Kasner equilibrium set 

In this subsection we perform a linearisation of the evolution system about the Kasner equilibrium points 



that form the parabola JC, given by Eqs. ( J106| ) and (112). We thus linearise Eqs. (113) 
the values 

E_=oS_(a;), = E^ = N x = E x = N_ = Q 



(120) about 



(129) 



where oS_(a;) is an arbitrary real- valued function of x. We have to treat the evolution equation ( |119|) for 
v in a special manner, due to the fact that the term E\ l d x il/il is singular on JC. We first linearise Eq. 



(118) for Q and then use the solution of this equation to show that the singular term in Eq. (119) can be 



16 The physical interpretation of this requirement is that matter does not affect the dynamics near the cosmological initial 
singularity. According to conjecture BKL I this condition will be satisfied except for special classes of models. 
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neglected when linearising it. Rescaling the time variable so that 7V"o 
system of linear ODE in t, with x-dependent coefficient qS_: 



dtE, 1 = 


2E-L 1 


d t £_ = 


-|(2- 7 )oS_n 


d t N x = 


2A x -(a cC £-)£i 1 


9 t S x = 


- 2V3o£_ S x 


cW- = 


2(1 + a/3oE_)JV_ 


d t n = 


|(2-7)(i + «^i)n 


d t v = 


i [ 3 7 - 2 - 3 (2 - 7 ) oSl ] v + 3 ^—^ qS_ TV 

7 


of Eq. (130) is 



^i 1 =6(x)e 2 



1, we then obtain the following 

(130) 
(131) 
(132) 
(133) 
(134) 
(135) 

(136) 
(137) 



We can use a reparametrisation (pq ) of x to set b{x) = 1. The resulting general solution of the linear 
ODE system is then given by (including the zeroth-order contribution to S_) 



(138) 

(139) 
(140) 



£_ 


= "**(*) 


fl " 1 


l + \k*{x)_ 


N x 


= [o 2 ( a; ) + ^t5 x ^)]e 2t 


Sx 


= o 3 ( I )e 2fc ( I » i 


iV_ 


= a^x)^ 1 -^ 1 





3 / 

= a 5 (x)e2( 2 - 


7)[l+|fc 2 (a>)]t 



o 6 (a;)e3l^- 2 -C 2 -T) fca (-)]*- 2 



fc(x) 



,2( 



^ 7 [l + i/fc 2 (a:)] 



02(2) + ^9x^(x) f i - 



3(2- 7 )[l + ifc2 (x)L 

For convenience, and to agree with the notation of Rendall and Weaver p3[ , we write 

oS_(a;) = --Lfc(x) 



(141) 
(142) 

(143) 

(144) 



for the limiting value of £_ as t — > — oo. The solution of the linear ODE system suggests that an arc JCa of 
the Kasner equilibrium set K attracts neighbouring orbits (i.e., is a local attractor). This arc is defined by 



the requirement that in each exponential function in the solution (138) - ( |143| ) the independent variable 

oo. The 



t has a positive coefficient, so that the solution approaches the equilibrium point (129J) as t 

size of the arc depends on whether the model is polarised (i.e., diagonal) or not, and on the equation of 

state parameter 7 of the fluid, as shown in Tab. || We stress that this linear analysis does not prove that 



Class of models 



Attracting arc JCa 



Vacuum/polarised 

Fluid/polarised 

Vacuum, or fluid with 1 < 7 < 2/unpolarised 



all of K. 

(37-2) 1 / 2 



V3(2- 7 ) 1/2 



<oE_ <0 



" 73^- 



<0 



Table 3: Attracting arc JCa on the Kasner parabola JC 

JCa is a local attractor. 

Over the past 11 years a number of rigorous analyses of the past asymptotic behaviour of G2 cos- 
mologies have been given which enable us to make precise statements about JCa- Firstly, Isenberg and 
Moncrief Q have proved that every polarised Gowdy vacuum solution with spatial topology T 3 is past 
asymptotic to a Kasner solution, showing that the Kasner equilibrium set JC is the global past attractor 
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for this class of models. Secondly, Kichenassamy and Rcndall 34 used an analysis based on the Fuchsian 
algorithm to prove that a general family| 17 | of unpolarised Gowdy vacuum solutions with spatial topology 
T 3 is past asymptotic to the arc Ka, as given in Tab. ||. In a recent development, Rendall [Q has used 
the /3-normalised scale-invariant FOSH evolution system for Gowdy vacuum spacetimes to argue that the 
arc K-a is in fact a local past attractor in the unpolarised case, for models which satisfy the so-called "low 
velocity" condition < i>Gowdy < 1, where the Gowdy "velocity parameter" corresponds to 



^Gowdy 



v^E^+E 2 ) 1 / 2 , (145) 



and thus quantifies the magnitude of the transverse shear rate of the timelike reference congruence erj. It 
is known, however, that the arc JCa is not the global past attractor for Gowdy vacuum spacetimes since 
solutions which develop so-called spikes violate the inequality < t>Gowd y < 1 at those points at which a 
spike occurs (see Rendall and Weaver |43], Berger and Moncrief [Bj and Hern and Stewart pqj). Finally, 
Anguige S has proved that a general family of diagonal Gi cosmologies with a perfect fluid matter source 
is past asymptotic to the arc JCa as given in Tab. g. 

It also follows from Refs. J33|, |M and [g| that the solution (137) - (143) to the linear equations gives 



the correct past asymptotic form of a general class of solutions in a neighbourhood of the local attractor 
JCa for polarised and unpolarised Gowdy vacuum spacetimes and for diagonal G2 cosmologies with a 
perfect fluid matter source. We anticipate that it will also do so for orthogonally transitive perfect fluid 
G2 cosmologies, but this remains to be proven. 

Finally, we note that if the restriction fE- > y=j- — -^tj^, which arises in the polarised perfect fluid 

G2 case in Tab. a, does not hold, then the peculiar velocity v of the fluid will remain significant as 
t — > — 00, hinting at the existence of another local attractor distinct from ICa- Experience with SH 
cosmologies (see, e.g., Ref. p0[) suggests that v will approach its extreme values, i.e., hmt^-ooV = ± 1. 
This matter requires further investigation. 

5.3 Asymptotic silence 

We now give a brief discussion of the conjecture BKL II, in the light of the previous two subsections. 
This conjecture is part of the folklore of mathematical cosmology and does not have a precise statement. 
We can best explain the essence of the conjecture by quoting from BKL B, p656: 

"... in the asymptotic vicinity of the singular point the Einstein equations are effectively 
reduced to a system of ordinary differential equations with respect to time: the spatial deriva- 
tives enter these equations 'passively' without influencing the character of the solution." 

Another way of expressing the idea heuristically is to say that the evolution at different spatial points 
decouples near the cosmological initial singularity. The FOSH format of the evolution system that we 
have given, namely 

A9 i X + B£ 1 1 3 I X = F(X), (146) 

(using the timelike area gauge with A/o = 1), sheds light on this idea of spatial decoupling, since we have 
shown that lim t _ > _ oc E\ l — when linit^-oo Q, = 0. This result means that as one follows a timeline into 
the past towards the cosmological initial singularity, the local null cone (and hence the local fluid sound 
cone embedded therein) collapses onto the timeline, showing that geometrical information propagation 
between neighbouring timelines is asymptotically eliminated, as illustrated in Fig. |l|. We shall refer to 
this phenomenon as "asymptotic silence" of the gravitational field dynamics as the cosmological initial 
singularity is approached. As regards the quotation from BKL ||, there are two ways of reducing the 



evolution system of PDE to a system of ODE. Firstly, one can set Ei 1 = in Eq. (146), obtaining the 
system of ODE that describes the dynamics on the unphysical boundary:]^] 

A9jX = F(X). (147) 

Secondly, one can consider the system of linear ODE in the neighbourhood of the Kasncr equilibrium set, 



Eqs. (130) - (136). We have seen that the system of linear ODE do produce the correct past asymptotic 

17 That is, a family whose initial data depends on four arbitrary real-valu ed f unctions. 

18 On the original idea of "silent cosmological models" , see Matarrese et al [B9| ; on their dynamical consequences , Ref. |20| . 

This is related to the so-called "velocity-dominated" system obtained by dropping the spatial derivatives from the 

evolution system but not from the constraint equations; see, e.g., Andersson and Rendall 
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Figure 1: Schematic representation of the phenomenon of "asymptotic silence" of the gravitational field 
dynamics in the approach of the cosmological initial singularity. 

form of the solutions near the locally attracting Kasner arc K-a- One can also consider the relation 
between Eqs. (147) and (146). In view of the fact that limt^__oo E\ l — 0, one might expect that the 
spatial derivative term BBi 9 X X in Eq. (|146| ) would be negligible compare d to the oth er two terms, as 
t — v — oo. Calculating each term in Eq. ( |146| ) using the asymptotic solution (137) - fll43| ) shows that the 



spatial derivative terms are in fact negligible asymptotically in a neighbourhood of the Kasner arc K,a- 
This property does not hold at isolated points in solutions which develop spikes, since certain partial 
derivatives in <9 T X become unbounded, so that the term B^ 1 d x X. is not negligible there. Nevertheless, 
the spatial derivative terms appear to act "passively", so that one still has silent, Kasncr-like dynamics 
locally. 



5.4 Isotropic initial singularities 



Our discussion in subsection 5.2 concerns the past asymptotic behaviour of general classes of G2 cosmolo- 
gies, which, according to conjecture BKL I, satisfy limt_ > _ 00 il = 0. There are, however, special classes of 
models that violate this condition, the most important being models with an isotropic initial singularity 
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(see Goode and Wainwright 23 ), which satisfy 



lim fi = 1 , lim v = , lim S' = 0. (148) 

t — > — oo t — > — oo t — > — oo 

Their evolution near the cosmological initial singularity is approximated by the flat FL model. We can 
understand the generality of the isotropic initial singularity by linearising the evolution equations ( [11 3D 
- ( |120| ) about the flat FL equilibrium point which lies on the unphysical boundary and is given by 

£l = (fl = l, = ^i 1 = S_ = N X = E X =N- =v . (149) 

Setting A/"o = 1, we thus obtain the following system of ODE: 

dtEx 1 = \{Z 1 -2)E 1 1 (150) 

9 t E_ = -§(2- 7 )E_ d t N x = i(3 7 -2)iV x (151) 

9 t E x = -| (2 - 7 ) E x d t N_ = ±(3 7 -2)7V_ (152) 

d t Vt = §(2- 7 )(l-fi) ft w = |(37 - 2) «. (153) 

Using a reparametrisation ( pq) of x, we obtain E\ x = exp(^(37 — 2)£), and 





= ai(a;)e^i (2 - 7)t 


7V X = ^(^e^ 37 " 2 )* 


(154) 


X 


= a 3 (a;)e-| (2 - 7)t 


7V_ = a 4 (x)e^ 3l - 2)t 


(155) 


n = 


= 1 + a 5 (x) e -i (2 - 7) * 


v = a 6 (x)e^- 2)t . 


(156) 



Note that it follows from the present result that for a G*2 cosmology to have an isotropic initial singularity 
the conditions 

= ai(x) = a^(x) = as(x) (157) 

need to be satisfied. This amounts to setting precisely half the initial data compared to the full G2 case. 
For rigorous results see Claudel and Newman |ll[ and Anguige and Tod ||. 

6 Concluding remarks and outlook 

In this paper we have shown how to formulate the EFE for orthogonally transitive G2 cosmologies with a 
perfect fluid matter source as an autonomous system of PDE with evolution equations in FOSH format, 
using scale-invariant dependent variables. As stated in the introduction, one of our goals is to provide 
a flexible framework for analysing Gi dynamics. A potential user of this paper, someone who wishes 
to apply the equation systems that we have derived to do numerical or rigorous analyses, need not be 
familiar with the orthonormal frame formalism. For such a reader the heart of the paper is the scale- 



invariant equation system in section |2.2| , the discussion of the gauge choices in section |3j, and the overview 
in section 



4.1. The explicit equation system in section 4.4 for the timelike area gauge may also be of use. 
On the other hand, for someone interested in the structure of the space of cosmological solutions of the 
EFE, the relevant sections are 4.2, 4.2 and a 



An important aspect of our formulation is that it incorporates the SH models as a special case, 
thereby shedding light on how SH dynamics influences G2 dynamics. The H- normalised scale- invariant 
dependent variables, when used in a dynamical systems setting, have proved effective in all four aspects 
of analysis of SH models, i.e., exact solutions, heuristic, numerical and rigorous. In particular, the 
dynamical systems framework suggested various heuristic ways of gaining insight into the SH dynamics, 
using local stability arguments and the notion of shadowing of orbits in the dynamical state space, in 
conjunction with the hierarchy of Bianchi invariant submanifolds. These methods have in turn led to 
proofs of various conjectures concerning SH dynamics (e.g., Ringstrom |44[| ) and to the prediction of 
new dynamical phenomena (e.g., asymptotic self-similarity breaking and Weyl curvature dominance; see 
Wainwright et al |35| ). Our initial success in describing some aspects of the G2 dynamics at early times in 
terms of a local past attractor, as given in section @, suggests that our new formulation of the evolution 
for perfect fluid G2 cosmologies will prove equally effective. The three main problems concerning G2 
dynamics that we intend to investigate in the future are 
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(i) the asymptotic dynamics at early times when the peculiar velocity is dynamically significant, 

(ii) the local stability of self-similar models and the asymptotic dynamics at late times, and 

(iii) the dynamics of G 2 cosmologies that are close to FL in some epoch. 

In conclusion, we believe that many of the ideas discussed in the present paper, when appropriately 
modified, will be of relevance in a much broader context in mathematical cosmology. For example, 
the notion of an infinite-dimensional scale-invariant dynamical state space with a hierarchical skeleton 
structure will be a useful guide for exploring more general cosmological spacetimes. We anticipate that 
concepts such as geometrical information propagation, asymptotic silence, and a past attractor located 
on the unphysical boundary of the dynamical state space, whose dynamics is described by a system of 
ODE, will play an important role in clarifying the dynamical content of the conjectures BKL I and BKL 
II concerning cosmological initial singularities. Likewise, these ingredients should be useful for studying 
almost-FL dynamical states near the cosmological initial singularity or at intermediate and late times, 
as well as other, more generic, aspects of these dynamical regimes. 
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A Appendix 

A.l Connection components in terms of frame variables 

The inverse area density of the G 2 -orbits is given by 

A' 1 = e 2 2 e 3 3 - e 2 3 e 3 2 . (158) 

Then it follows from the dimensional commutator equations ( |l2| ) - (|j]) that 

)u-i dtei1 ■ i dxN nrcm 

a = -N — u\ = ei — — - (159) 

ei 1 N 



r _ 1 d t A ! ! d x A 

~A Sl = ^ ei T 



= \N^ U -^ ai = ~W^T (160) 



<*- = -^AN- 1 (e 3 3 d t e 2 2 -e 2 2 d t e 3 3 + e 2 3 d t e 3 2 -e 3 2 d t e 2 3 ) (161) 

" x = ^/ 3 ^ l 1 (e 3 3 d x e 2 2 -e2 2 d x e 3 3 + e 2 3 d x e 3 2 -e 3 2 d x e2 3 ) (162) 

CT * = ^ 3 -AN- 1 (e 2 3 d t e 2 2 -e 2 2 d t e 2 3 + e 3 2 d t e 3 3 -e 3 3 d t e 3 2 ) (163) 

n - = ^^ e i l (e2 3 d x e 2 2 -e 2 2 d x e 2 3 + e 3 2 d x e 3 3 -e 3 3 d x e 3 2 ) (164) 

fil = \AN^{e 2 3 d t e 2 2 -e 2 2 d t e 2 3 -e 3 2 d t e 3 3 + e 3 3 d t e 3 2 ) (165) 

n+ = -±Aei 1 (e 2 3 d x e 2 2 -e 2 2 d x e 2 3 -e 3 2 d x e 3 3 + e 3 3 d x e 3 2 ). (166) 

A. 2 Scale-invariant curvature variables 

We define additional /3-normalised curvature variables by 

(0*, S..., £..., «... ) := (-|*R, *S..., £..., H... )/(3/3 2 ) . (167) 
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Then we obtain for orthogonally transitive Gi cosmologies with a perfect fluid matter source the following 
expressions 



N on- zero 3-Ricci curvature variables: 



n k = -*(E 1 1 d IB -r-%A)A+(N$+Nl) 

S+ = -liE^dx-^A+KNl+Nl) 

S- = \{E 1 1 d x -r-2A)N x +lN+N_ 

S x = - | (E^ d x -r-2A)N_ + l N+ N x . 



(168) 
(169) 
(170) 
(171) 



Non-zero characteristic Weyl curvature variables f~\ 



(£-±n x . 



- I (E 1 1 d, - r) A + § (TV 2 + iV 2 ) + ± £+ - ± (S 2 . + £ 2 ) + ± 7G; 1 ilv 2 (172) 
-7V_£_-7V X I] X (173) 

± I (Eh 1 d x - r + 3S+ - A) (£_ ± JV X ) T f iV+ (E x T M) + § (2_ - A7V X ) (174) 



(fxTtt-) = ±i(^i 1 ^- r T3S+-^)(S x TiV_)±|iV + (E_±Ar x ) + i(i] x+ ^iv_)(175) 

A. 3 Line element and scale-invariant dependent variables for area gauges 

Introducing in the separable area gauge a line element of the form 

.2 



ds 2 =ll[- e 2 **-*) di 2 + e 2 ^*^ dx 2 + e 2Mot m(x) (e p(t ^ (dy + Q(t, x)dz) 2 + e - p ^ dz 2 
the area expansion rate is given by 

Then we obtain the following expressions for the non-zero scale-invariant dependent variables: 

(AT 1 , Si 1 ) = (A/- -\A/- -V-s) 

(£+, U) = (l(l-Aro 1 d t g),E l 1 d m f) 

(E_,JV X ) = _^(^-ia iP) - El x d x P) 

(E X ,7V_) = 1 c p (AA - 1 a f Q, Si^Q) 



A = - 



2\/3 

1 dlnm(a;) t 



d.i: 



^l 1 , 



(176) 
(177) 

(178) 
(179) 
(180) 
(181) 

(182) 



and R = — v3 S x and N + = V3 iV_. Employing an x-reparametrisation (pq) to set m(x) = exp(— 2D SSi x), 
with Z? sa a constant, implies A — D sa E\ 1 . Apart from the sign of t, these expressions reduce to the Gowdy 
vacuum line element of Rendall and Weaver [Q when A/o = |, m(x) = 1, /(£, x) = | [A(t, x) + 3i] and 
g (t, x) = \ [ \(t , x) - t ] . Note that then d x E-J^= 0. 

A. 4 Propagation of constraint equations 

Propagation of dimensional constraint equations: 

N' 1 d t (C com ) A 12 = -(a + /3 + V3a_)(C com ) A i2-(V3a x +r! 1 )(C com ) A 3i 



N- l d t {c com y 
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-{a + p + V3a_) (C com ) A i2 - {V3a 

— ^2 (Ccodacci)l 

-{a + 13- V3a_) (C com ) A 3 i - (V3a 

— S3 (Ccodacci)l 

N- x d t {C G&u6 ) = -(a + / S)(C G aufi)-4(Mi + ai)(Ccodacci)i 

A''~ 1 9t(Ccodacci)l = ~ (tt + 3/3) (Ccodacci)l - 3 (Ml - Ol) (CGaufi) 

M- l d t {C p ) = - (1 - 3E+) (Cfc) . 



fii)(C, 



com) 12 



(183) 

(184) 
(185) 
(186) 
(187) 



3 Again, we correct some sign errors in the expressions given in Refs. WN and |18|, 
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For simplicity, in deriving these results we have used the gauge fixing condition (|12J ) as an identity. 
Propagation of dimensionless gauge fixing conditions: 

N^dtiCjjU = ( 9 + l)(C [> ) nc (188) 

Mo 1 dtiCfiU = ( 9 + 3S+)(C c> ) sa (189) 

AA- 1 a t (C & ) fc = [(g + 3S+)+3( 7 -l)(l-S + )](C £> ) fc (190) 

M- 1 d t {C ) sync = ( g + 3S+)(C c> ) sync . (191) 
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